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Abstract 

We analyze the potential of the CERN Large Hadron Collider (LHC) to study anomalous quartic 
vector-boson interactions ZZ77, ZZZ^, W + W~jj, and W + W~ Z7 through the weak boson 
fusion processes qq — > ^577 and qq — > qqjZ(-^ with t = e or \i. After a careful study of the 

backgrounds and how to extract them from the data, we show that the process pp — > is 
potentially the most sensitive to deviations from the Standard Model, improving the sensitivity to 
anomalous couplings by up to a factor 10 4 (10 2 ) with respect to the present direct (indirect) limits. 
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I. INTRODUCTION 



Within the framework of the Standard Model (SM), the structure of the trilinear and 
quartic vector-boson couplings is completely determined by the SU{2) L x U(l)y gauge 
symmetry. The study of these interactions can either lead to an additional confirmation of 
the model or give some hint on the existence of new phenomena at a higher scale The 
triple gauge-boson couplings have been probed at the Tevatron 2] and LEP 3], 0] through 
the production of vector-boson pairs, however, we have only started to study directly the 
quartic gauge-boson couplings If any deviation from the SM predictions 

is observed, independent tests of the triple and quartic gauge-boson couplings can give 
important information on the type of New Physics (NP) responsible for the deviations. For 
example, the exchange of heavy bosons can generate a tree level contribution to four gauge- 
boson couplings while its effect in the triple-gauge vertex would only appear at one one-loop, 
and consequently be suppressed with respect to the quartic one. Further information on the 
NP dynamics can also be provided by determining whether NP reveals itself in the form 
of anomalous four-gauge couplings involving only weak gauge bosons or in those involving 
photons or in both. 



ous^ couplings arises 
Due to phase 



At present the scarce experimental information on quartic anoma 
from the processes e + e~ — > W + W~ r y, Zjj, ZZ'j, and z/z/77 at LEP 
space limitations, the best sensitivity is attainable for couplings involving photons which 
should appear in the final state. Photonic quartic anomalous couplings can also affect 77Z 
and 'yyW productions at Tevatron j^, |t| , however, it was shown in Ref . |s| that even with an 
integrated luminosity of 2 fb _1 , the Tevatron experiments can only probe the gauge quartic 
couplings at the level of precision obtained at LEP. In the near future, both photonic and 
non-photonic quartic gauge couplings will be tested in pair production of gauge bosons at the 
LHC via weak boson fusion (WBF) jslQ]- In the long term, high sensitivity to anomalous 
photonic fo_ couplings is expected at a next e+e" Unear co„ ld e r flQ, as we,, as 
high energy 77 |L2J, LLj] , and [la] colliders. 

In this work, we study the potential of the LHC to probe the photonic quartic vertices 
ZZ77, W + W~ , yy, W + W~ Zj, and ZZZ'j. The motivation for this study is two folded: first, 
even at LHC energies, the best experimental sensitivity is expected for couplings involving 
photons due to phase space limitations. Second, if a signal is observed, the comparison of 
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the processes here studied, which are only sensitive to photonic quartic operators, with the 
observations for processes also dependent on non-photonic couplings, such as weak gauge 
boson pair production, could reveal some symmetries of the underlying dynamics. 

We perform a detailed analysis of the most sensitive channels that are the production via 
WBF of photon pairs accompanied by jets, i.e., 

p + p->q + q->j+3+l + l, (1) 

and the WBF production of a pair of jets plus a photon accompanied by a lepton pair, where 
the fermions originate from the decay of a Z° or a virtual photon i. e. 

p + p -»• qq -»• j + j + 7 + (Z* or 7* ->) t + f , (2) 

with i = e or //. The advantage of WBF, where the scattered final-state quarks receive 
significant transverse momentum and are observed in the detector as far-forward/backward 
jets, is the strong reduction of QCD backgrounds due to the kinematical configuration of 
the colored part of the event. 

The process depicted in Eq. (J2J) receives contributions from all four-gauge-boson vertices 
that we are interested in, while only the ZZ77 and W + W~^^ vertices are relevant for the 
process in Eq. ((TJ). We previously studied the reaction (JTJ) in Ref. Q|. Here, we revisit the 
limits there obtained after taking careful account of the QCD uncertainties in the background 
evaluation and analyzing strategies to minimize it, and compare them with the expected 
sensitivity from (J2J). 

This paper is organized as follows. We present in Sect. II the effective operators we 
analyzed in this work. Section III contains our analysis of the signal and backgrounds, as 
well as the attainable limits ate the LHC We draw our conclusions in Sect. VI. 



II. EFFECTIVE QUARTIC INTERACTIONS 

We parameterize in a model independent form the possible deviations of the SM pre- 
dictions for the photonic quartic gauge couplings with the assumptions that NP respects 
SU(2)l x U(l)y gauge invariance and that no new heavy resonance has been observed. 
In this scenario the gauge symmetry is realized nonlinearly by using the chiral Lagrangian 
approach as in Ref. p] . Following the notation of Ref. ^| , the building block of the chiral 
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Lagrangian is the dimensionless unimodular matrix field T,(x) 

.Lf a (x)T a 



cxp 



(3) 



where v = {V2G F y l . The fields are the would-be Goldstone fields and r a (a = 1, 2, 3) 
are the Pauli matrices. The SU(2)l x U(1)y covariant derivative of £ is defined as 

r3 



(4) 



We focused our attention on genuine photonic quartic interactions, i.e. the new inter- 
actions do not exhibit a triple gauge boson vertex associated to them. In our framework, 
genuine quartic operators appear at next-to-leading order (0(p 4 )), however, there is no 
genuine photonic quartic interaction at this order. Therefore, we considered the next order 
(0(p 6 )). There are 14 effective photonic operators which respect SU{2) C custodial symmetry 
as well as C and V, 

,2 



9' 
A 2 



A 2 " 



99 
A 7 



^Tr(W^W I ')Tr(y Q! K) + ^Tr(^^ a )Tr(^V;) + k^Ti(W^V a )i:i{W^V a ) 



+k™Tr(W^V u )Tr(W» a V a ) + k^Ti(W^V a )Ti(W^V v ) 
k h Q Ti{B^B^)Ti{V a V a ) + h h c Tt{B^B^)Tz{V v V a ) 



+ 



+ 



+k b 1 Tr(B llu V a )Tr(B^V a ) + k h 2 Ti{B^V v )Ti{B^V c 
k^Tr{W^B» v )Tr{V a V a ) + k^TiiW^B^TiiV^) + k^^W^V^^V, 
+k^Tr(W liU V v )Tr(B^ a V a ) + kfT^{W lii y a )Tr{B^V v 



(5) 



where V» = £t, B^ = r 3 5^/2, and = T a W^/2, with B„ v and W% being 

respectively the U(l)y and SU(2)l field strength tensors. Here, e is the electromagnetic 
coupling, g = e/ sin 6*^ = e/s w , and g' = g/c w with c w = \J\ — s^. A is a mass scale 
characterizing the NP. 

It is interesting to express the effective interactions in (jSJ) in terms of independent Lorentz 
structures. The lowest order effective W + W~ r y r y and interactions are described in 

terms of four Lorentz invariant structures 



Wn 7 



eV 



F^F^W +a W- , 



F^F^{W +U W- + W~ U W+) 



(6) 
(7) 
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ZS = -^F^Z a Z a , (8) 
Z2 = - e ^F^F» a Z»Z a , (9) 



while the lowest order effective the ZZZ^ interactions are given by 



Z o = -^-F^Z a Z a , (10) 
Z c = -^-F^Z^Z»Z a . (11) 



The remaining W + W Z7 interactions are parameterized as 

W Z = -e 2 g 2 F^Z^W +a W~ , (12) 
2 2 

Wf = -^-F^Z^iW+'W- + W~ V W+) , (13) 
2 2 

Wf = -^Z-F^w+W-Z* + W~W+Z a ) , (14) 

2 2 

W 2 Z = -^ J -F^{W+ M W- a Z„ + W- a W +a ^) , (15) 
zc w s w 

2 2 

W 3 Z = -^-F^(W+ a W-Z" + W- a W u + Z a ) . (16) 

The Fevnman rules for the quartic couplings induced by the above operators can be found 
in Ref. 

Eq. (J3J) can be conveniently rewritten in terms of the above independent Lorentz struc- 
tures, neglecting possible 4W, 4Z, WW ZZ as well as Goldstone boson vertices, as 

h 1 hi h 1 h 1 

£= f^ + w2) + f 2 {z2 + w2) + f 2 z^ + ^z2 (17) 

h z k z k w 

with 

k] = k™ + k\ + k? for % = 0, c, 1 , (18) 

^23 = ^2 ^2 ^2™ ^3 + k™ , (19) 

^0 = — (*o + *a ) - —(^0 + *i) + c^(C + fcD , (20) 

c w 

k z c = ^(k™ + + fc») - — (k b c + k b 2 ) + c zw (k™ + K + k™) , (21) 

s w c w 

UW _ 2^_h w - —k b 4- r C99^ 
"'O — "'O ' ^zw^o j 
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kV = k? + h™ for i = 1,2,3 



(23) 



(24) 



and c zw — (c^, — s^) / (2c w s w ) . 

Before we study the phenomenological consequences of anomalous quartic vertices, we 
should stress that the effective Lagrangian (fT7|) can also be obtained using a linear repre- 
sentation of the SU{2) L x U(l)y gauge symmetry with the presence of a Higgs boson in 
the spectrum However, in this case, the lowest order terms that can be written are 
of dimension 8 and they lead to different relations between the couplings associated to the 
independent Lorentz structures. Moreover, they generate both photonic and non-photonic 
genuine quartic vertices whose strength is in general related, unlike in the non-linear case. 

III. SIGNALS AND BACKGROUNDS 

In this work we study the reactions (JIJ and (J2J) at the LHC We evaluated numerically 
the helicity amplitudes of all the SM subprocesses leading to the jj'77 and jj"fl + l~~ final 
states where j can be either a gluon, a quark or an anti-quark in our partonic Monte Carlo. 
The SM amplitudes were generated using Madgraph |lq in the framework of Helas [17| 
routines. The anomalous interactions arising from the Lagrangian © were implemented as 
subroutines and were included accordingly. We consistently took into account the effect of all 
interferences between the anomalous and the SM amplitudes and did not use the narrow- 
width approximation for the vector boson propagators. We considered a center-of-mass 
energy of 14 TeV and an integrated luminosity of 100 fb -1 for LHC. 

It is important to note that the operators in Eq. © lead to tree-level unitarity violation 
in 2 — > 2 processes at high energies 0. The standard procedure to avoid this unphysical 
behavior of the cross section and to obtain meaningful limits is to multiply the anomalous 
couplings (kj) by a form factor 



where m 77 is the invariant mass of the final state photon pair in subprocesses like ZZ — > 77 
and WW 77. For subprocesses of the type ZZ — > Zj — > / + Z~7 and WW — * Z7 — > 




(25) 
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the anomalous couplings are multiplied by a form factor 

(2 \ ~ n 

1 + xH, (26) 

where m^- 7 is the invariant mass of the final state lepton pair plus a photon. Of course 
using this procedure the limits become dependent on the exponent n and the scale A u which 
is not longer factorizable. In our calculations, we conservatively choose n = 5 and A u = 2.5 
TeV for the LHC. 

At e + e~ colliders the center-of-mass energy is fixed and the introduction of the form 
factors (j23|) and (J2lj|) is basically equivalent to a rescaling of the anomalous couplings k\, 
therefore we should perform this rescaling when comparing results obtained at hadron and 
e + e~ colliders. For example, the LEP limits should be weakened by a factor ~ 1.6 for our 
choice of n and A u . 

Altogether the cross sections for processes (|T} and (J2J) can be written as 

k j - k j2 

O = Osm + (Jintcr + -^J- O"ano , (27) 

where <r sm , o" int er, and cr ano are, respectively, the SM cross section, interference between the 
SM and the anomalous contribution and the pure anomalous cross section. 

A. p + p^-j + j + 7 + 7 

This process receives contributions from ZZ77 and WWj'y vertices which get modified 
by all operators in Eq. (JSJ). However, as seen in the first line in Eq. (fTTj) there are only four 
independent Lorentz invariant structures contributing to this process which, consequently, 
is only able to give information on the four linear combinations of anomalous couplings 
corresponding to the four coefficients, k] , {i = 0, c, 1, 23) defined in Eqs. (fTHj) and (fTTJj) . 

Process receives contributions from W* and Z* productions in association to photons 
as well as from WW and ZZ fusion processes 

p + p^q + q + (W* + W* or Z* + Z*) ^g + g + 7 + 7. (28) 

In order to reduce the enormous QCD background we must exploit the characteristics of the 
WBF reactions. The main feature of WBF processes is a pair of very far forward/backward 
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tagging jets with significant transverse momentum and large invariant mass between them. 
Therefore, we required that the jets should comply with 

p^ 2) > 40 (20) GeV , |tfe (1 , a) | < 5.0 , 

\Vh ~Vj 2 \ > 4 - 4 , Vh-Vh<Q and (29) 
ARjj > 0.7 . 

Furthermore, the photons are central, typically being between the tagging jets. So, we 
require that the photons satisfy 

£^ (1 ' 2) >25GeV , |r/ 7(12) | < 2.5 , 

min^,^} + 0.7 < r? 7(1 2) < max^, r/ J2 } - 0.7 , (30) 
Ai2j 7 > 0.7 and Ai? 77 > 0.4 . 

Further reduction of the SM background can be achieved by a cut in the invariant mass 
distribution of the 77 pairs. As illustrated in Fig. the invariant mass distribution for 
the SM background contribution is a decreasing function of the 77 invariant mass while the 
anomalous contribution first increases with the 77 invariant mass reaching its maximum 
value at m 77 ~ 1000 GeV and then decreases. Consequently, in order to enhance the WBF 
signal for the anomalous couplings we imposed the following additional cut in the diphoton 
invariant mass spectrum 

400 GeV < m 77 < 2500 GeV. (31) 

We present in Table U] the values for cr ano for each of the independent linear combination 
of anomalous couplings in Eqs. (fTHj) and (fl~9"j) after applying the cuts in Eqs. (j2~9"j) - (j3"Tj) . 

These results were obtained using \/I as the factorization scale in the parton distribution 
functions. We have further assumed a 85% detection efficiency of isolated photons, leptons 
and jet-tagging. With this the efficiency for reconstructing the final state j + j ' + 7 + 7 is 
(0.85) 4 ~ 52% which is included in the results presented in Tables HI and HT1 The interference 
terms (cri nter ) between the anomalous and SM amplitudes turn out to be negligible. As 
expected the WW fusion process due to Wq (WJ) leads to a larger anomalous contribution 
(by a factor ~ 2.5) than the ZZ fusion ones due to the Zq (Z]). 

The evaluation of the SM background (cr sm ) deserves some special care since it has a 
large contribution from QCD subprocesses whose size depends on the choice of the renor- 
malization scale used in the evaluation of the QCD coupling constant, a s (/i_R), as well as 
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FIG. 1: Normalized invariant mass distribution of the 77 pair for the reaction pp — ► 77JJ. 



Coupling Constant 


OWo (pb x GeV 4 ) 




2.1 x 10 7 


kj 


1.5 x 10 6 


kj 


6.0 x 10 6 


h 1 

ft 23 


4.3 x 10 5 



TABLE I: Results for <r ano for process Eq. (J2) [see Eq. (|27j)] . We considered n = 5 and A u = 2.5 
TeV; see Eq. (|25|). All results include the effect of the cuts in Eq. (|29[). (|3T))) . and (|3*T|) as well as 
photon detection and jet-tagging efficiencies. 

on the factorization scale fip used for the parton distribution functions. To estimate the 
uncertainty associated with these choices, we have computed cr sm for two sets of renormal- 
ization scales, which we label as /im,2(0> an d for several values of //p. is defined 

such that Q!g(/iRi(£)) = a s (^p^ 1 )a s (^pt^ 2 ) where and are the transverse momentum of 
the tagging jets and £ is a free parameter varied between 0.1 and 10. The second choice of 
renormalization scale set is //R2(£) = with being the subprocess center-of-mass 

energy. 

In Table ITT1 we list <r sm for the two sets of renormalization scales and for three values of the 
factorization scale jjp = V^, \/I/10, and p^ in where p^ in = min(^,p^). As shown in this 
table, we find that the predicted SM background can change by a factor of ~ 8 depending on 



9 





0"sm (ft>) 




HR = /WO 


Mil = /WO 




[IF = V$ 


HF = Pmin 


[if = vi/10 


[IF = V§ 


= Pmin 


Hf = vi/10 


0.10 


3.2 


5.3 


4.1 


1.3 


2.2 


1.7 


0.25 


2.2 


3.6 


2.8 


1.1 


1.9 


1.4 


1.00 


1.4 


2.4 


1.9 


0.91 


1.5 


1.2 


4.00 


1.1 


1.8 


1.4 


0.78 


1.3 


1.0 


10.0 


0.94 


1.6 


1.2 


0.71 


1.2 


0.96 



TABLE II: Results for cr sm for process Eq. see Eq. (|27j) and text for details. All results 
include the effect of the cuts in Eq. (|2*9*|) . (|3l7|) and (pTTj) as well as photon detection and jet-tagging 
efficiencies. 

the choice of the QCD scales. These results indicate that to obtain meaningful information 
about the presence of anomalous couplings one cannot rely on the theoretical evaluation of 
the background. Instead one should attempt to extract the value of the SM background 
from data in a region of phase space where no signal is expected and then extrapolate to 
the signal region. 

In looking for the optimum region of phase space to perform this extrapolation, one 
must search for kinematic distributions for which (i) the shape of the distribution is as 
independent as possible of the choice of QCD parameters. Furthermore , sin ce the electroweak 
and QCD contributions to the SM backgrounds are of the same order |l8j, this requires that 
(ii) the shape of both electroweak and QCD contributions are similar. Several kinematic 
distributions verify condition (i), for example, the azimuthal angle separation of the two 
tagging jets which was proposed in Ref. [l9| to reduce the perturbative QCD uncertainties 
of the SM background estimation for invisible Higgs searches at LHC. However, the totally 
different shape of the electroweak background in the present case, renders this distribution 
useless. 

We found that the best sensitivity is obtained by using the 77 invariant mass. As can 
be seen in Fig. [TJ the shape of the SM distribution is quite independent of the choice of the 
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0.1 



0.08 - 



0.06 



0.04 - 



0.02 








FIG. 2: Ratio R(£) defined in Eq. (|32*|) for the process pp — > jjjj at LHC. 

QCD parameters. As a consequence most of the QCD uncertainties cancel out in the ratio 

cx(400 GeV < m 77 < 2500 GeV ) 



R(0 



(32) 



cx(100 GeV < m 77 < 400 GeV ) 
This fact is illustrated in Fig.|2]where we plot the value of the ratio R(£) for different values of 
the renormalization and factorization scales. The ratio R is almost invariant under changes 
of the renormalization scale, showing a maximum variation of the order of ±6% for a fixed 
value of the factorization scale. On the other hand, the uncertainty on the factorization 
scale leads to a maximum variation of 12% in the background estimation. We have also 
verified that different choices for the structure functions do not affect these results. 

Thus the strategy here proposed is simple: the experiments should measure the number of 
events in the 77 invariant mass window 100 < m 77 < 400 GeV and extrapolate the results 
for higher invariant masses using perturbative QCD. According to the results described 
above we can conservatively assign a maximum "QCD" uncertainty (QCD unc ) of ± 15% to 
this extrapolation. 

In order to estimate the attainable sensitivity to the anomalous couplings we assume 
that the observed number of events is compatible with the expectations for = 1) and 

\xp = so the observed number of events in the signal region coincides with the estimated 
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number of background events obtained from the extrapolation of the observed number of 
events in the region where no signal is expected; for this choice the number of expected 
background events is iV back = <r sm £ where C stands for the integrated luminosity. For an 
integrated luminosity of 100 fb^ 1 for LHC, this corresponds to iVback = 143. Moreover, we 
have added in quadrature the statistical error and the QCD uncertainty associated with the 
backgrounds. Therefore, the 95% limits on the quartic couplings can be obtained from the 
condition 



k 



2 



iV a no = ^ x C x cr ano < lMy/N^ + (iV back x QCD unc ) 2 . (33) 

For the sake of completeness we show the results on the expected sensitivity using purely 
statistical errors and for two values of QCD unc : our most conservative estimate [15 %), 
and a possible reduced uncertainty (7.5 %), which could be attainable provided NLO QCD 
calculations are available. Assuming that only one operator is different from zero, so no 
cancellations are possible, we find 

\k™' b ' m /A 2 \ < 3.3 (3.9) [4.8] x lO" 6 GeV~ 2 , 
\k™' b ' m /A 2 \ < 1.3 (1.5) [1.8] x 10- 5 GeV~ 2 , 

\k™' b ' m /A 2 \ < 6.2 (7.2) [8.9] x lO" 6 GeV" 2 , (34) 
\k™' b ' m /A 2 \ < 2.3 (2.7) [3.3] x lO" 5 GeV~ 2 , 
\k^' m /A 2 \ < 2.3 (2.7) [3.3] x 10~ 5 GeV" 2 . 

We notice that the constraint on k^ b,m and k^' m are exactly the same as they both modify 
in the same form and amount the process (0) as seen in Eq. ()17|) . 

Finally, let us comment that the limits on fc^ 1 ' 6 '" 1 / A 2 and k™' b ' m / A 2 can be directly trans- 



lated on constraints on the coefficients ao, c , of the operators introduced in Ref. 
substitution a 0c = 4g 2 kQ C (see Eq. (IT%|0 . 



12] with the 



b. p + P -> j + j + 7 + 1+ + r 

This process receives contribution from the four-gauge coupling vertices ZZZj and 
WWZ'y as well as from ZZ r y r y and WW r ) r ). We have imposed a minimal set of cuts to 
guarantee that the photons, charged leptons and jets are detected and isolated from each 
other: 

pj (2) > 40(20) GeV , p £ T > 25 GeV , Ej, > 25 GeV , 
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\r ]l/ \<2.5 , \r] j(12) \ < 5.0 , 

\Vh-Vj 2 \> 4A > Vh-Vj2<0, (35) 
mi HVji,Vj 2 } + 0- 7 < < max{?7 il , r) j2 } - 0.7 , 
^Rjj ijijt) > 0.7 , ARg+e-^g-j > 0.4 . 

Furthermore, in order to single out the events containing Z° bosons and to enhance the 
WBF signal for the anomalous couplings ZZZ'y and WW Z^ we have imposed the following 
additional cuts on the lepton-lepton (m«) and lepton-lepton-photon (m 7 «) invariant masses: 

\m u -M z \<20 GeV and 400 GeV < m jU < 2500 GeV. (36) 

In Table lllll we display the values of <7 ano after cuts for each anomalous couplings k\ in 
Eq. (jHJ), with \ip = y/l. These results include the effect of detection and tagging efficiencies; 
85% efficiency for detecting isolated photons and leptons and for tagging a jets. With this, 
the efficiency for reconstructing the final state j ' + j ' + 7 + l + l~ is (0.85) 5 ~ 44%. We have 
added the contributions from final states containing electrons and muons. Once again, we 
verified that the interference terms cr inter are neglig ible. 

A detailed study of the results in terms of the different Lorentz structures involved show 
that the invariant mass cut on the lepton-lepton invariant mass suppresses the contribu- 
tions from the W + W~ r y r y Lorentz structures Wq and W] in relation to those containing the 
VVZ'-f and quartic vertices (V = W or Z). However, we find that none of the Lorentz 

structures involving these vertices is clearly dominant and that there are important inter- 
ference effects between the different Lorentz structures contributing to the same anomalous 
operator, which are the order of 10-30% and can be destructive or constructive. 

The evaluation of the SM background in this case is also subject to QCD uncertainties, 
as in the previous reaction. We found that for our reference value /im(£ = 1) an d = V« 

<j sm = 0.10 fix (37) 

Changes in the factorization and renormalization scales, can modify this prediction by 
a factor ~ 5. Thus, again, the best strategy to accurately determine the sensitivity to the 
anomalous coupling is to extract the value of the SM background from data in a region of 
phase space where no signal is expected and then extrapolate to the signal region. Following 
the discussion in the previous section, we find that the £ + £~7 invariant mass distribution 
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Coupling Constant 


^ano (pb x GeV 4 ) 


ft/Q 


4.6 x 10 7 


b w 
fv c 


9.2 x 10 6 


kf 


2.9 x 10 7 


h w 
«, 2 


1.3 x 10 7 


kf 


1.0 x 10 7 


k b 


6.9 x 10 6 


k\ 


1.9 x 10 6 


k\ 


4.7 x 10 6 


k\ 


1.8 x 10 6 





1.1 x 10 7 


kf 


3.2 x 10 6 




9.0 x 10 6 




4.3 x 10 6 




3.6 x 10 6 



TABLE III: Results for 

°"ano f° r the process (J2j); see Eq. (|27j) . cr ano is obtained for the anomalous 
coupling kj/A 2 in units of GeV -2 . We considered n = 5 and A u = 2.5 TeV; see Eq. (|26j) . 

is suitable to estimate the SM background and reduce the QCD uncertainties. We have 
defined the ratio 

(x(400 GeV < m« 7 < 2500 GeV ) 
^ ~ o (100 GeV < m« 7 < 400 GeV ) ' ( ' 

and evaluated the behavior of R(£) under changes of the renormalization and factorization 
scales. We determined that R(£) can be known within an accuracy of ± 15% when we use 
leading order calculations. 

In order to extract the attainable limits on the anomalous couplings we assumed a lu- 
minosity of C = 100 fb -1 and that the observed number of events is compatible with the 
expectations for /iRi(£ = 1) and //p = y/S, i.e. the expected number of background events 
in the signal region is iVb ac k = 10. We have added to the statistical error associated with 
this background the theoretical error associated to the uncertainty in the extrapolation of 
the background. However, given the limited statistics, the sensitivity is dominated by the 



14 



statistical error. The 95% CL. constraints on the anomalous couplings are 
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which have been obtained including a 15% QCD uncertainty, However, to the precision 
quoted, the impact of this uncertainty is minimal. 

Comparing the limits in Eqs. ()39j) with the corresponding ones from the process ([TJ 
in Eq. we see that, despite the limited statistics, the presence of the VVZj vertex 
(V = W or Z) makes the process pp — > jjjl + l~ most sensitive to the presence of NP leading 
to anomalous four-vector couplings which respect the SU(2)i x Z7(l)y gauge invariance as 
well as the SU(2) C custodial symmetry. One of the reasons for the process pp jj r yl + l~ 
to be more sensitive to anomalous interactions is that almost all Lorentz structures lead 
similar contributions and that more Lorentz structures contribute to this reaction than in 
pp — > jj'77 for a given effective operator. 

One must keep in mind, however, that the results in Eqs. (}3"4*|) and (}3*9"j) were obtained 
under the assumption that only one operator is different from zero, so no cancellations were 
possible. If cancellations are allowed the process (JIJ may become the most sensitive one to 
the presence of the relevant photonic quartic operators. 
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IV. DISCUSSION 



We are just beginning to test the SM predictions for the quartic vector boson interactions. 
Due to the limited available center-of-mass energy, the first couplings to be studied should 
contain photons. In particular, the direct searches at LEPII have lead to constraints of 
the order |^| < 0(1CT 2 GeV~ 2 ) for the couplings in Eq. (J3J), and no significantly better 
sensitivity is expected from searches at Tevatron. Anomalous quartic couplings contribute 
at the one-loop level to the Z physics via oblique corrections as they modify the W, 
Z, and photon two-point functions. Consequently they can be indirectly constrained by 
precision electroweak data to < 0(1O -4 GeV~ 2 ). 

Higher energy colliders will be able to test quartic gauge couplings involving photons 
as well as to probe non-photonic vertices WW' (V,V = W or Z) Even at LHC 

energies, due to phase space limitations, the best experimental sensitivity is expected for 
couplings involving photons which can be part of the final state. Moreover, in the event that 
a departure from the SM predictions is observed, inference of the underlying dynamics can 
only be obtained by comparing the observations in different channels, for instance between 
those involving triple and quartic-gauge couplings. In this respect it will also be important to 
know whether NP reveals itself in the form of anomalous four-gauge couplings involving only 
weak gauge bosons or in those involving photons or in both. For instance, in the framework 
of chiral lagrangians, where no light Higgs state is observed, the photonic four vertices are 
expected to be suppressed with respect to the non-photonic ones since they appear one order 
higher in the momentum expansion. An anomalous signal only in the photonic couplings 
could indicate that there are additional symmetries forbidding the non-photonic vertices. 

With this motivation, in this work we have analyzed the production of two jets in asso- 
ciation with a photon pair, or with a photon and a £ + £~ pair at LHC as tests of anomalous 
bosonic quartic couplings involving one or two photons. In this study we have taken careful 
account of the theoretical uncertainties associated with the evaluation of the SM background. 
We have proposed the best strategy to estimate the expected SM background by extrapola- 
tion of the data taken in a region of phase space where no signal is expected, minimizing the 
theoretical uncertainty associated to this extrapolation. The final sensitivity to the different 
couplings is given in Eqs. (JMj) and (|39j). In particular, we found that in the framework of 
SU(2)l x U(l)y gauge invariant NP in which the deviations from the SM prediction for 
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the W77 vertices are related to the strength of the anomalous VVZ'y vertex, the process 
PP ~^ 33ll + l~ is the most sensitive to all possible operators, despite the limited statistics, 
barring possible cancellations. It can lead to constraints \^\ ^5 1.2-6.3 x 10~ 6 GeV -2 . 

In conclusion we have shown that the study of the processes (0) and (j2J) at LHC can 
test quartic anomalous couplings that are four orders of magnitude weaker than the existing 
limits from direct searches and two orders of magnitude weaker than any indirect constraints. 
It is interesting to notice that if no signal is found the LHC will lead to limits that are similar 
to the ones that could be attainable at an e + e~ collider operating at yfs = 500 GeV with a 
luminosity of 300 fb" 1 2, Q- 
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